Precise QCD predictions on top quark pair production 
mediated by massive color-octet vector boson at hadron colliders 
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We present a theoretical framework for systematically calculating next-to- leading order (NLO) 
QCD effects to various experimental observables in models with massive COVB in a model in- 
dependent way at hadron colliders. Specifically, we show the numerical results for the NLO QCD 
corrections to total cross sections, invariant mass distribution and Afb of top quark pairs production 
mediated by a massive COVB in both the fixed scale (top quark mass) scheme and the dynamical 
scale (top pair invariant mass) scheme. Our results show that the NLO QCD calculations in the 
dynamical scale scheme is more reasonable than the fixed scheme and the naive estimate of the NLO 
effects by simple rescaling of the LO results with the SM NLO K-factor is not appropriate. 



In many extensions of the Standard Model (SM), mas- 
sive Color-Octet Vector Boson (COVB) is necessarily en- 
gaged at the TcV scale, for example, in the top-color 
warped (RS) or universal extra dimensions @, 13] , techni- 
color 0| and chiral color models In all these cases, 
the COVB could have large impacts on the interaction 
of top quarks, which are being copiously produced at 
the CERN Large Hadron Collider (LHC). With a large 
sample of ti data, the CDF collaboration at the Teva- 
tron has recently reported an observation of a large 
Forward-Backward asymmetry (Afb) in tt production, 
= 0.158 ± 0.075, compared with the SM prediction 
0.058 ± 0.009 The disagreement is more profound 

in the region of large ti invariant mass, where CDF re- 
ported A FB (m tf > 450GcV) = 0.475±0.114, and the SM 
gives 0.088 ± 0.013 @. This leads to a more than 3.4 a 
deviation from the SM prediction @. Similarily, the DO 
collaboration [l(| has also reported the total Afb to be 
Af* = 0.196 ± 0.065 using 5.4fb _1 of data. Further- 
more, D0 has also measured the charge asymmetry of 
the leptons from top decay, A l FB = 0.152 ± 0.04. These 
results have generated extensive theoretical studies on 
this observable in various models beyond the SM. Among 
these, models with massive COVB are in particular at- 
tractive, c./., Ref. [ll[ and references therein. There have 
also been substantial efforts in searching for the signal of 
COVB at the Tevatron and LHC, which can shows up as 
a cle ar reso nant peak in the tt invariant mass distribu- 
Whilc current experimental limits depend 
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tion ^ 

on the detailed choices of couplings EMI, they never- 
theless indicate that COVB with mass below 1 TeV is 
severely constrained. 

It is well known that QCD effects play an important 
role in ti production. The NLO QCD corrections to SM tt 
production, which significantly enhance the ti total cross 
sections, have been calculated for a long time (mm 



In the SM, Afb is related to higher order QCD radia- 
tion effects, which first appear at 0{a 3 s ) 0. Complete 
NLO corrections to this observable are not available cur- 



rently, but calculation based on soft gluon resummation 
indicates that higher order QCD effects are small 22 1. 
Recently electroweak corrections to Afb have also been 
calculated and are found to slightly increase the asym- 
metry 23[. In the case of massive COVB, its QCD gauge 
interaction is uniquely determined by its color content, 
resembling a SM gluon. Therefore, it's reasonable to ex- 
pect that higher order QCD effects will also have signifi- 
cant impacts on processes mediated by massive COVB, at 
least at an energy scale comparable to the mass of COVB. 
This has motivated the model dependent calculation of 
COVB production by gluon fusion [24[ and the model in- 
dependent (using dimension-six operators) calculation of 
ti production mediated by COVB J2||. However, a com- 
plete NLO analysis of the QCD effects to models with 
massive COVB in the resonant ti region is still absent. 
In this letter, we present the model independent com- 
plete NLO QCD corrections to top quark pair produc- 
tion mediated by a general massive COVB, and show 
the detailed numerical analysis of top quark pair produc- 
tion, including invariant mass distribution and Afb at 
the NLO level. We also show that the NLO corrections 
significantly stabilize the renormalization and factoriza- 
tion scale dependence, as compared to the LO results. 

Below, we briefly outline our approach to systemati- 
cally calculating the NLO QCD effects to processes of 
COVB production. We consider a model independent 
massive COVB originated from a broken SU(3) gauge 
group. The effective Lagrangian for color-octet vector 
G M in unitary gauge can be written as 
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Tr G fiV G' J ' u + MqTt G^G^, 



(1) 



where a = 1,..,8 are the broken SU(3) "color" indices, 
and fi — 0, ..,3 are the Lorentz indices. G^ = G a ^JT a — 



1 While this work was completed, Ref. |26| appeared, which calcu- 
lates the singly production of color-octet vector boson in narrow 
width approximation. 
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{d^Gl - d v G a ^T a is the field strength tensor, where T a 
is the conventional Gellmann matrix with the normaliza- 
tion Tr[T a T b ] = \5 ab . The mass of the COVB is denoted 
by Mq. The QCD color interaction between COVB and 
SM gluon can be easily implemented by changing the 
ordinary derivative into covariant derivative: 



dfj,G 



g 1 f abc A b G c in 



(2) 



where g\ is the coupling constant of QCD. The La- 
grangian in Eq. (JTJ), after the replacement in Eq. @, is 
already invariant under the conventional SU(3) C trans- 
formation. If desired, NLO QCD calculation can be 
done with the Feynman rule derived from the above La- 
grangian, where unitary gauge is chosen for the broken 
SU(3) gauge symmetry. However, it is well known that 
loop calculation in unitary gauge is inconvenient because 
of the violent ultra-violet (UV) behavior of the propa- 
gator. Instead, we choose to carry out the calculation 
in conventional 't Hooft-Feynman gauge for the broken 
gauge group. To this end, we separate the longitudinal 
component (the would-be Goldstone boson) of the COVB 
from Eq. ((T|), by modifying the mass term in Eq. ((T|) as 
follows, 



G^x) -> G^x) = U{xy I -dp + G^(x)j U(x), (3) 

where we have introduced the 7r field (the would-be Gold- 
stone field) via U(x) = e I7r I?, and the symmetry 
breaking scale / = Mclgn, with gi being the coupling 
constant of the broken SU(3) gauge symmetry. It is easy 
to check that the mass term is now invariant under the 
gauge transformation: 

G^x) V(x) (g^x) + j-d^ V( X y, (4) 

U(x) -> V(x)U(x), (5) 

where V(x) is a finite gauge transformation V(x) = 
exp(ia2(x)T a ). Similar to unitary gauge, the interaction 
between G M , 7r and QCD gluon is obtained by changing 
the ordinary derivative into covariant derivative, as in 
Eq. ([2]). Using Eqs. ([2]) and ([3]), the classical Lagrangian 
for COVB turns into 



-^Tr G^G^ 



(6) 



d u ^D u 



Note that the classical Lagrangian £q is non- 
renormalizable, i.e., when expanding U(x), it contains 
operators of dimension larger than 4, whose coefficients 
are suppressed by at least one power of /. Since we 
are only interested in low energy QCD effects of COVB 
rather than its UV-completion, these higher dimensional 
operators can be safely neglected. Expanding the mass 
term in Eq. we find that there is a kinetic mixing 



between the COVB and the would-be Goldstone boson: 
Mc'K a {x)d tl G'^ l {x). This mixing can be canceled by in- 
troducing a gauge fixing term: 

F 2 = ~{d»Gl{x) + M G * a {x)) 2 , (7) 

which is similar to the gauge fixing term for the QCD 
Lagrangian in 't Hooft-Feynman gauge: 



Fi = -^a;( x )) 2 . 



The corresponding ghost Lagrangian is given by 

SF 

Cg = ^ {x) Jeix) u1 ^ 



(8) 



(9) 



where i = 1 for the QCD ghost, and i — 2 for the broken 
SU(3) ghost. 6i(x) are the infinitesimal gauge transfor- 
mation parameters for the corresponding gauge group. 
The complete Lagrangian for the gauge sector now can 
be written as 



C = Cr 



-QCD 



F 1 +F 2 +C q , 



(10) 



where Cqcd is the classical QCD Lagrangian. All the 
Feynman rules determined by gauge symmetry can then 
be derived from Eq. (fTU|) . and couplings between COVB 
and SM fcrmions can be introduced through vector or 
axial-vector current interaction. We have checked that 
this set of Feynman rules are in agreement with those 
derived in Ref. [53] for the case of RS model • 

Before continuing, we should define the precise mean- 
ing of the NLO QCD corrections in this paper. At LO, 
there are three parts of contributions to the tt cross sec- 
tion in models with massive COVB: the squared SM am- 
plitudes, the interference between the new physics (NP) 
amplitudes and the SM amplitudes, and the squared NP 
amplitudes. The NLO QCD corrections in this paper 
refer to the 0(a s ) corrections to these three parts, re- 
spectively. All the relevant Feynman rules can be de- 
rived from the effective Lagrangian given in Eq. (jlOp . 
The results obtained in this way reflect the model inde- 
pendent corrections from QCD interaction. As a final 
comment, we note that our NLO QCD corrections are 
the non-Abelian analogy of the QED corrections to 
boson production at hadron collision J5ij], and the ne- 
glected contributions in our calculation are similar to the 
genuine weak corrections there [5j| , which can be mostly 
absorbed into redefinition of boson and fcrmion couplings 
at LO. 

Extending the approach shown in our paper [53], we 
calculate the one-loop renormalized helicity amplitudes, 
with the unstable particle treated in the complex mass 
scheme [5^] . The loop integrals with complex arguments 
appearing in the one-loop amplitudes are evaluated with 
ONELOOP [30J1 . The real emission matrix elements are 



generated by a modified version of MADGRAPH [31 
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Soft and (or) collincar singularities of the real correc- 
tions are dealt with by the dipole method [32[ imple- 
mented in the MADDIPOLE package [33[. Throughout 
our calculation, the pole mass of top quark is chosen as 
tth = 173.1 GeV, and a s (M z ) = 0.118. LO cross sec- 
tions are obtained with the CTEQ6L parton distribution 
function (PDF) [3~ij with one- loop running of a s , while 
NLO cross sections are obtained with the CTEQ6M PDF 
with two- loop running of a s . In the numerical calculation 
below we present the results for a benchmark axial-gluon 
model [35[ resulting from a simutaneous fit of ti total 
cross section, Afb, and dijet resonance searches. In this 
model, The coupling between the massive COVB and 
quarks are chosen as 



v q {m t ) 
vt{m t ) = 



= 0, 
0, 



a q {m t ) 
at{m t ) = 



2.2, 
-3.2, 



(11) 



where giv q ( t ) and g\a q (t) are the vector and axial cou- 
pling of the light quark (top quark) to the COVB, de- 
fined at the scale mt- The evolution for f q,t (a*) un- 
der the change of scale is given by gi{y)v{fj) q ^ = 

gMv(vo) q Aas(ri/aM) 15/i2M M, where /3 = 
23/3 is the first QCD beta- function coefficient for N c = 3 
and rif = 5, and the evolution equation for a (Jj t(/Lt) has 
the same form. Unless specified, the mass of the COVB 
is chosen as 2 TeV. We note that while such choice of 
parameters is at the margin of currently allowed param- 
eter space 36l~l38|. the qualitative behavior of the NLO 
QCD effects arc quite general and can be applied in more 
realistic models. For the calculation of accurate numeri- 
cal results for other parameter choices, the corresponding 
results can be directly obtained from our Fortran code. 

We define the NP cross section, <tnp, as the difference 
between ti cross section in a model with a massive COVB 
and the SM: 



cr NP = OsM+NP — 0SM, 



(12) 



where the SM cross sections, including both the qq- and 
gg-channels, are calculated with the program MCFM Q. 

In Fig. [TJ we plot the NP contribution to the cross 
section at the LHC with s/s = 7 TeV, as a function 
of Mq- The bands reflect the scale uncertainties esti- 
mated by varying the renormalization (fi r ) and factoriza- 
tion (fJ.f) scales around their default values by a factor 
of 2 with [i r = fif. We present the results in two bench- 
mark schemes, namely the fixed scale scheme (FSS), i.e. 
the scales are fixed to be m t , and the dynamical scale 
scheme (DSS), i.e. the scales are set to be the invari- 
ant mass of top quark pair m t t- We find that the NLO 
QCD effects are large and small in the FSS and the DSS, 
respectively, for our choice of parameters and the naive 
estimate of the NLO effects by simple rescaling of the 
LO results with the SM NLO K-factor is not appropri- 
ate. It is also clear that the inclusion of NLO QCD effects 
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FIG. 1. onp, defined in Eq. (|12[> . as functions of Mg for 
FSS (left) and DSS (right). The black (red) bands are the 
LO (NLO) uncertainties, estimated by varying the scales 
around their default values by a factor of two within each 
scheme. The blue dashed lines are the naive estimates of the 
NLO effects by simple rescaling of the LO results with the 
SM K-factor. The couplings are chosen as in Eq. 



strongly reduces the theoretical uncertainty of the NLO 
cross section in the both scheme, compared with the LO 
ones. 
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FIG. 2. LO (dashed lines) and NLO (solid lines) invariant 
mass distribution of the top quark pair at the LHC with ^/s = 
7 TeV are shown in the upper panel. In the bottom panel 
the differences between two different schemes at the LO and 
NLO are shown. In both panels, red curves and black curves 
correspond to FSS and DSS, respectively. The couplings are 
chosen as in Eq. 



Fig. [5] gives the LO and NLO invariant mass distribu- 
tions of the top quark pair at the LHC with ^fs = 7 TeV, 
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including contribution from SM LO and NLO top quark 
pair production. It can be seen that the NLO results in 
the DSS has relatively small corrections compared with 
the LO one, while the NLO QCD corrections in the FSS 
significantly change the shape of the LO curve, leading 
to the reduction of the width of the resonance. This is 
due to the fact that the NLO width can be expressed 
analytically as 



A FB by 3 - 4%. 



NLO 



(M) 



40/167 15 M£\ 

TT \ 12 4 (JL 2 J 



(13) 

From Eq. (JT3J), it is obvious that the width of COVB 
is significantly reduced at NLO in FSS, while the large 
logarithmic contributions are absent in DSS. The pre- 
dictions in the two schemes at NLO are close to each 
other, while at LO they show large difference, as shown 
in the lower panel of Fig. [2l The difference between the 
two schemes reflects the uncertainties of the theoretical 
prediction. Hence, the NLO result leads to a smaller 
theoretical uncertainties in m t i distribution, which could 
improve the accuracy of extracting the theory parame- 
ters of NP models from comparing to experimental data. 
We also note that similar conclusion holds in other cases, 
e.g., KK gluon in RS model. 

From the calculations of the total cross section and the 
invariant mass distribution, we see that the perturbation 
expansion may break down when FSS is chosen because 
the COVB resonance mass is usually about an order of 
magnitude larger than the top quark mass, leading to 
very large logarithms in the perturbation series. In con- 
trast, the scale in DSS is related to the dynamics of the 
process and thus is a more reasonable choice. 

Fig. El shows the LO and NLO contributions to the 
Apb as a function of the mass of the COVB at the Teva- 
tron with y/s = 1.96 TeV in the center of mass frame of 
tt pair. The results are given for both the total asym- 
metry (bottom bands) and the asymmetry in the large 
invariant mass region (upper bands), m t j > 450 GeV, 
respectively. The bands in Fig. [3] reflect the scale uncer- 
tainties of the theoretical predictions, which are obtained 
by simutaneously varying (i r and (if in the numerator 
and denominator of A-pB around m t by a factor of 2 with 
(i r = (if = (i. Here we choose the FSS, because the SM 
results Apg = 0.158 ± 0.075 is also obtained at the scale 
(X = rn t . In both cases, the most significant effects of 
the NLO corrections are the reduction of scale depen- 
dences. It is worth pointing out that the numerator of 
Afb is dominated by NP contributions, while the denom- 
inator is dominated by SM tt cross section. Therefore at 
LO the large scale dependence in the numerator, which is 
mainly from scale dependence of NP coupling v q j {(J.) and 
a Qi t((i), can not be canceled by the corresponding scale 
dependence in the denominator. The predictions for Afb 
at (i = m t arc depicted in Fig. [3] as dashed lines. It can 
be seen that the NP contributions at NLO only reduce 
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FIG. 3. NP contributions to the Afb in the tt center of mass 
frame as a function of the mass of the COVB at the Tevatron 
with y/s = 1.96 TeV. In the upper lines only the Afb in the 
large invariant mass region, mti > 450 GeV, are plotted, while 
in the bottom total Afb are shown. The bands reflect scale 
dependence, estimated by varying the scales around their cen- 
tral values by a factor of 2. The predicitions for (i — mt 
are depicted in dashed lines. The couplings are chosen as in 
Eq. (Till. 

Figs. [1] [2l and [3] show the cross sections as functions 
of Mq, invariant mass distribution and the Afb as a 
function of Mq for a COVB with axial couplings, re- 
spectively. For the Afb, it would also be interesting to 
illurstrate the results for a COVB with both vectorial 
and axial couplings. For this purpose, we consider an- 
other set of couplings for a COVB with Mq = 3 TeV as 
follows. 



v q {m t ) 
vt(m t ) 



-1, 
-3, 



a q {m t ) 
a t {m t ) - 



(14) 



The results for this set of couplings are shown in Fig. @] 
Again we can see that the NP contributions at NLO 
only reduce Afb by a small amount for default scale 
choice {(i = m t ). 

As a further application of our results, we plot in Fig. [5] 
the NP cross section, ctnp, for a specific RS model consid- 
ered in search for resonance decaying into top quark pairs 
by ATLAS collaboration [39[, where the left-handed and 
right-handed couplings (modulo g s ) between KK gluon 
and light quarks are chosen as 

c L , q (M G ) = c R>q {M G ) = -0.2, (15) 

and similarily the couplings between KK gluon and top 
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FIG. 5. LO and NLO predictions for <jnp in a specific RS 
model [3!j . The expected and observed limits on cross section 
are extracted from the Ref . . 



FIG. 4. NP contributions to the Afb in the ti center of mass 
frame as a function of the mass of the COVB at the Tevatron 
with -y/s = 1.96 TeV. In the upper lines only the Afb in the 
large invariant mass region, m t j > 450 GeV, are plotted, while 
in the bottom total Afb are shown. The bands reflect scale 
dependence, estimated by varying the scales around their cen- 
tral values by a factor of 2. The predicitions for n — mt 
are depicted in dashed lines. The couplings are chosen as in 
Eq. (Ell). 

quark are chosen as 

C£,t(M G ) = l, c R , t {M G ) = A. (16) 

Also plotted in Fig. [5] are the experimental exclusion 
limit extracted from the same paper. Our exact theoret- 
ical predictions arc given for three different scale choices, 
fx = mt, fx = Mq and fi = ra t i, at both LO and NLO. 
It can be seen that the NLO QCD effects are moder- 
ate in the fi = m t i case, while they are large in both the 
fx = mt and fx = Mq cases. We also plot in Fig.[5]the LO 
results (the black solid line) with only the contributions 
from NP squared amplitudes, i.e., without the interfer- 
ence with the SM amplitudes. It is clear that considering 
the NP squared amplitudes contributions alone obviously 
underestimates the NP cross sections, and including the 
interference contributions is necessary for reliable predic- 
tions. We note that while the mass limit for KK gluon is 
very different at LO for the three kinds of scale choices, 
their differences arc significantly reduced at NLO. Hence, 
the NLO results can be used for precise extraction of mass 
limit for KK gluon. 

Finally, we give the numerical results of the NP cross 
section for the same set of parameters as in Fig. with 
Mq fixed to be 1.5 TeV. The results are presented for 
three different scale choices, i.e. fi = m t , /i = Mq and 
fx = m t i with [i r = fif = /i. The LO (NLO) results 



are calculated with CTEQ6L (CTEQ6M) in Table U and 
with MSTW2008LO (MSTW2008NLO) in Table HI re- 
spectively In the LO calculations, we only estimate the 
theoretical uncertainties by varying the renormalization 
and factorization scales around their central values by a 
factor of 2. In the NLO calculations, we also consider 
the PDF uncertainties, and the first error is scale uncer- 
tainty and the second is PDF uncertainty. It can be seen 
from Tables U and |TT] that scale uncertainties are reduced 
from LO to NLO in all three scale choices, and the dif- 
ferences between different scale choices are also reduced 
when going from LO to NLO. Moreover, the predictions 
of total cross sections for CTEQ and MSTW PDFs are 
obviously different at LO, but consistent within the PDF 
uncertainties at NLO. 
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TABLE I. The NP cross sections for top quark pair produc- 
tion mediated by KK gluon at the LHC (y/s = 7 TeV) at 
the LO (CTEQ6L) and NLO (CTEQ6M). Theoretical un- 
certainties are estimated by varying the renormalization and 
factorization scales around their central values by a factor of 
2. 

In conclusion, we have presented a theoretical frame- 
work for systematically calculating NLO QCD effects to 
various experimental observables in models with massive 
COVB in a model independent way. Specifically, we show 
the numerical results for the NLO QCD corrections to to- 
tal cross sections, invariant mass distribution and Afb of 
top quark pairs produced by mediating a massive COVB. 
Our results show that, for our choice of parameters, the 
NLO corrections of the NP cross section are much larger 
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TABLE II. The NP cross sections for top quark pair produc- 
tion mediated by KK gluon at the LHC (^/s — 7 TeV) at the 
LO (MSTW2008LO) and NLO (MSTW2008NLO). Theoreti- 
cal uncertainties are estimated by varying the renormalization 
and factorization scales around their central values by a factor 
of 2. 



in the FSS than in the DSS, but the NLO QCD calcu- 
lations in the dynamical scale scheme is more reasonable 
than the fixed scheme, and the naive estimate of the NLO 
effects by simple rescaling of the LO results with the SM 
NLO K-factor is not appropriate. Moreover, for invari- 
ant mass distribution, we find that the NLO QCD cor- 
rections change the width of the resonant particle in the 
DSS slightly, and the differences in the results using FSS 
and DSS are reduced after including NLO QCD effects. 
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